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(Dated: February 7, 2008) 
We report branching fraction measurements for exclusive decays of charged and neutral 
B mesons into two-body final states containing a charmonium meson. We use a sample of 
22.72 ± 0.36 milUon BB events collected between October 1999 and October 2000 with the 
BABAR detector at the PEP-II storage rings at the Stanford Linear Accelerator Center. The 
charmonium mesons considered here are J/^ , ijj{2S), and Xci, and the light meson in the decay 
is either a K, K* , or tt". 

PACS numbers: ll.30.Er, 13.25.Hw 



I. INTRODUCTION 

Decays of B mesons to two-body final states containing 
a charmonium resonance (J/f/', iPi'^S), Xci) constitute 
a very sensitive laboratory for the study of electroweak 
transitions, as well as the dynamics of strong interactions 
in heavy meson systems. In particular, neutral B decays 
to these final states are expected to exhibit a significant 
CP asymmetry, the magnitude of which is cleanly related 
to standard model parameters iQ. 

The tree level and leading penguin diagrams for the 
decay modes we consider are shown in Fig. |l|. Due to 
the contributions of non-perturbative QCD interactions 
in the final state, assumptions must be made in estimat- 
ing the expected branching fractions of these modes, and 
therefore these estimates have some degree of model de- 
pendence. A number of such estimates hav e ap peared in 
the literature |, | |, | |, @, |, |, 0, g. The one 
model-independent element common to all of these pre- 
dictions is the requirement from isospin symmetry that 
the ratio of the charged to neutral partial widths should 
be unity, and that this should hold separately for each 
light meson accompanying the charmonium meson in the 
final state. 

Here we report the measurement of branching fractions 
of B mesons to a charmonium resonance accompanied by 
a kaon or tt*' meson. The channels measured are listed 
in Table H. Here and throughout this paper for each final 
state mentioned its charged conjugate is also implied. We 
reconstruct J/t/j decays to lepton pairs £^£^, where £ is 
either an electron or muon. 

Our large data sample permits a measurement of these 
branching fractions with a precision superior to previous 
experiments. The simultaneous measurement of a num- 
ber of final states allows us to determine ratios such as 
vector to pseudoscalar kaon and heavy to light charmo- 
nium states production. Many systematic errors cancel 
when these ratios are extracted from a single data set us- 
ing very similar event selection criteria, further increasing 
the usefulness of our results for the validation or devel- 
opment of phenomenological models. 

Another highly relevant input for the understanding of 
strong interactions in B decays is the measurement of po- 
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1: Leading Feynman diagrams for the decays we con- 



larization in vector-vector final states, which is reported 
in another publication |l3|. Finally, the branching frac- 
tion of B —f J/tp 7r+ is measured using a specific analysis 
method, reported in [Q . 



II. THE BABAR DETECTOR 

The BABAR detector is located at the PEP-II e+e~ 
storage rings operating at the Stanford Linear Acceler- 
ator Center. At PEP-II, 9.0 GeV electrons coUide with 
3.1 GeV positrons to produce a center-of-mass energy of 
10.58 GeV, the mass of the T{4:S) resonance. 

The BABAR detector is described elsewhere ||l^; here 
we give only a brief overview. Surrounding the interac- 
tion point is a 5-layer double-sided silicon vertex tracker 
(SVT) which gives precision spatial information for all 
charged particles, and also measures their energy loss 
(dE/dx). The SVT is the primary detection device for 
low momentum charged particles. Outside the SVT, a 
40-layer drift chamber (DCH) provides measurements of 
the transverse momenta px of charged particles with re- 
spect to the beam direction. The resolution of the pt 
measurement for tracks with momenta above 1 GeV/c is 
parameterized as: 



*Also with Universita di Perugia, Perugia, Italy. 
tAlso with Universita della BasiUcata, Potenza, Italy. 



Pt 



= 0.13pt( GeV/c)% + 0.45%. 



(1) 



TABLE I: Branching fractions and decay modes considered in 
this paper. We always reconstruct the J/xj) in the t^l~ decay 
mode. 



Branching fraction 




Secondary decay 


measured 




modes used 


B" ^ J/^PK" 


K'' ^ 


K^h K'^ ^ 7T+7T- or n'^ Tv'* 




K'' -. 


Kl 


B+ -^ J/-ipK+ 


~ 




B'' -> J/i> K"-' 


^*o_ 


* K+ TT- or K^s tt"; K^s -^ tt+tt- 


B+ -^ J/ipK'-+ 


K*+ - 


^ K+ tt" or Kl 7V+; K^ ^ ^+^- 


B'^^J/^n" 


- 




B" -^ i^i2S)K'i 


V(2S) 


-*£+£- or J/ij TT+TT"; 

Kl ^ Tv+n- 


B+ -^ ^{2S)K+ 


^(2S) 


-> ^+^- or i,{2S) -^ J/4> TT+TT- 


B -^ Xci^s 


Xci -^ 


J/V- 7; Kl -^ TT+TT^ 


B+ ^ x^iK+ 


Xcl -» 


J/V- 7 


B" - XciK"' 


Xcl -> 


J/t/; 7; 7^*" ^ K+ TV' 



The drift chamber also measures dE/dx with a precision 
of 7.5%. Beyondtheouter radius of the DCH is a detector 
of internaUy reflected Cherenkov radiation (DIRC) which 
is used primarily for charged hadron identification. The 
detector consists of quartz bars in which Cherenkov light 
is produced as rclativistic charged particles traverse the 
material. The light is internally reflected along the length 
of the bar into a water-filled stand-off box mounted on 
the rear of the detector. The Cherenkov rings expand 
in the stand-off box and are measured with an array 
of photomultiplicr tubes mounted on its outer surface. 
A CsI(Tl) crystal electromagnetic calorimeter (EMC) is 
used to detect photons and neutral hadrons, as well as 
to identify electrons. The resolution of the calorimeter is 
parameterized as: 



oiEl 
E 



2.3% 



(^(GeV))4 



® 1.9%. 



(2) 



The EMC is surrounded by a superconducting solenoid 
that produces a 1.5-T magnetic field. The instrumented 
flux return (IFR) consists of multiple layers of resistive 
plate chambers (RPC) interleaved with the flux return 
iron. In addition to the planar RPC layers in the flux 
return, there is an additional cylindrical layer just outside 
of the EMC. The IFR is used in the identiflcation of 
muons and neutral hadrons. 

Data acquisition is triggered with a two-level system. 
The first level (Level 1) monitors trigger information from 
the DCH and EMC, and generates a trigger upon de- 
tection of track or cluster candidates. The second level 
(Level 3) retains events in which the track candidates 
point back to the beam interaction region (L3 DCH trig- 
ger), or EMC clusters candidates remain after the sup- 
pression of hits which have less energy than a minimum 
ionizing particle or arc uncorrclated in time with the rest 
of the event (L3 EMC trigger). Over 99.9% of BB events 
pass either the L3 DCH or L3 EMC trigger. A fraction of 
all events that pass the Level 1 trigger are passed through 



Level 3 to allow monitoring of the Level 3 trigger perfor- 
mance. 



III. DATA SAMPLE 

The data used in these analyses were collected between 
October 1999 and October 2000 and correspond to an 
integrated luminosity of 20.7 fb^^ taken on the T{AS) 
and 2.6 fb~ taken off- resonance at an energy 0.04 GcV 
lower than the peak, which is below the threshold for BB 
production. The data set contains 22.72 ± 0.36 million 
BB events. 



IV. COORDINATE SYSTEM AND REFERENCE 

FRAMES 

We use a right-handed coordinate system with the z 
axis along the electron beam direction and y axis up- 
wards, with origin at the nominal beam interaction point. 
Unless otherwise stated, kinematic quantities are calcu- 
lated in the rest frame of the detector. The other refer- 
ence frame we commonly use is the center of mass of the 
colliding electrons and positrons, which we will call the 
center-of-mass frame. 



V. PARTICLE RECONSTRUCTION 

The reconstruction of exclusive B decays begins with 
identifying candidates for the decay products. Charged 
particles are reconstructed as tracks in the SVT and/or 
DCH. Leptons and kaons are identified with information 
from the DCH, the EMC (for electrons) the IFR (for 
muons), and the DIRC (for kaons). Photons are identi- 
fied based on their energy deposition in the EMC, and 
K^ are identified from either energy deposition in the 
EMC or a shower in the IFR. 



A. Track Selection 

In general, tracks used in this analysis are required to 
include at least 12 DCH hits to ensure that their mo- 
menta and dE/dx are well measured. In addition, tracks 
are required to have pt > 100 MeV/c, and to point back 
to the nominal interaction point within 1.5 cm in xy and 
3 cm in z. Roughly 95% of the solid angle about the 
interaction point in the centcr-of-mass frame is covered 
by 12 or more DCH layers. 

We make exceptions to this requirement for two 
types of particles: pions from ii'!?, which do not origi- 
nate at the nominal interaction point, and pions from 
V'(25') -^ J/tp 7r+7r~, which frequently do not have suf- 
ficient transverse momenta to traverse 12 layers of the 
DCH. Any track found in the DCH or SVT is used in 
reconstructing these particles. 



TABLE II: Summary of electron identification criteria. Variables used are: AE /Ax, the energy loss measured in the DCH; 
E/p, the ratio of the EMC cluster energy to the momentum measured in the tracking spectrometer; Nays, the number of EMC 
crystals forming the cluster; LAT, the lateral energy distribution pM of the EMC cluster; A42, one of the Zernike moments 
[ [L7[ of the EMC cluster; and 6c, the Cherenkov angle measured in the DIRC. In addition, the fraction of electrons in inclusive 
J/ip events that pass each set of criteria is shown, along with the fraction of pions with momentum above 1 GeV/c that pass 
the selection requirements. 





DCH-only 


Loose 


Tight 


Very tight 


dE/dx (measured-expected) 


-2 to +4 (Jmcas 


-3 to +7 (Jmcas 


-3 to +7 (Tnicas 


-2 to +4 (Tnioas 


E/p 


- 


0.65 - 5.0 


0.75 - 1.3 


0.89 - 1.2 


Ncrys 


- 


>3 


>3 


>3 


LAT 


- 


- 


0.0 - 0.6 


0.1 - 0.6 


A42 


- 


- 


- 


<0.11 


9c (measured-expected) 


- 


- 


- 


-3 to +3 (Tnioas 


Efiiciency (%) 


94.9 


97.2 


95.4 


88.2 


TV misID (%) 


21.6 


4.8 


1.2 


0.1 



TABLE III: Summary of muon identification criteria. Variables used are: -Eemc, the energy deposited by the muon candidate 
in the EMC (this requirement is only applied for tracks within the fiducial coverage of the EMC); A^iayers, the number of IFR 
layers with hits; Nx, the number of nuclear interaction lengths traversed; \Nx — Nx{exp)\, the difference between the number of 
nuclear interaction lengths traversed and the expectation for a muon of the measured momentum; (A'hit), the average number 
of hits per IFR layer; RMShit, the RMS of the distribution of the number of hits on each layer; /hit, the fraction of layers 
between the innermost and outermost hit layers that also have hits (this requirement is only applied in the region covered 
partly or entirely by the endcap IFR system, 0.3 < 9 < 1.0); Xifri the x^ oi the track in the IFR; and Xmatchi the x^ of the 
match between the IFR track and the track from the central detector. In addition, the fraction of muons in inclusive J/tp events 
that pass each set of criteria is shown, along with the fraction of pions with momentum above 1 GeV/c that pass the selection 
requirements. 





MIP 


Very Loose 


Loose 


Tight 


Very tight 


^EMC ( GeV) 


< 0.5 


< 0.5 


<0.5 


0.05 - 0.4 


0.05-0.4 


Mayers 


- 


> 1 


> 1 


> 1 


> 1 


iVA 


- 


>2 


>2 


> 2.2 


> 2.2 


\Nx-Nx{exp)\ 


- 


< 2.5 


<2.0 


< 1 


<0.8 


(iVhit) 


- 


< 10 


< 10 


< 8 


< 8 


RMShit 


- 


<6 


<6 


< 4 


<4 


/hit 


- 


>0.1 


> 0.2 


> 0.3 


>0.34 


XlFR 


- 


- 


< 4 X A'^Iayors 


< 3 X Mayers 


<C X iViayers 


Xmatch 


- 


- 


< 7 X A^laycrs 


< 5 X Mayers 


< 5 X Mayers 


Efiiciency (%) 


99.6 


92.2 


86.2 


70.3 


67.0 


n misID (%) 


57.9 


14.5 


7.0 


2.4 


2.1 



B. EMC cluster reconstruction 



Photon Candidate Selection 



The energy deposited in contiguous crystals of the 
EMC is summed into a cluster. The distribution of en- 
ergy among the crystals is used to discriminate between 
clusters arising from electromagnetic and hadronic show- 
ers. The variables used to describe this distribution are 
the lateral energy (LAT) ||l^ and the Zernike moments 
Amn flj- LAT is a measure of the radial energy pro- 
file of the cluster; the Zernike moment A42 measures the 
asymmetry of the cluster about its maximum. Electro- 
magnetic showers have LAT peaked at about 0.25 and 
A42 close to zero, while showers from hadrons have a 
broader distribution in LAT, and extend to larger values 
of A42 . 



Photons are identified as EMC clusters that do not 
have a spatial match with a charged track, and that have 
a minimum energy of 30 MeV. To reject clusters arising 
from noise hits, LAT is required to be less than 0.8. 



D. Electron and Muon Identification 

We derive substantial background rejection from the 
positive identification of electrons and muons within the 
sample of charged tracks. For electrons, the variables 
that distinguish signal from background include LAT and 
A42, the ratio of energy measured in the EMC to mo- 
mentum measured in the tracking spectrometer (E/p), 
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dE/dx measured in the DCH, and the Cherenkov angle 
9c measured in the DIRC. 

For identifying muons, the presence of an energy depo- 
sition consistent with a minimum ionizing particle in the 
EMC, and the details of the distribution of hits in the 
IFR are used. In particular, the number of interaction 
lengths traversed in the IFR N\ must be consistent with 
expectations for a muon, both the average and variance 
of the number of hits per layer must be small, and the 
fit of a track to the hits must have low x^, both within 
the IFR (Xipi^) and in the match between the IFR and 
central detector track (Xmatch)- 

Since the optimal tradeoff between efficient selection 
and suppression of backgrounds varies between decay 
modes, there are several sets of criteria used to select 
leptons. These are defined in Table M for electrons 
and Table [II for muons. In addition to these criteria, 



we also restrict the lepton selection to a fiducial region 
within which the efficiency is well-known from control 
samples, and the material in the detector is accurately 
modelled in the Monte Carlo. The accepted range in po- 
lar angle is 0.410 < 6 < 2.409 r ad for electrons and 
0.30 < 6 < 2.70 rad for muons. This corresponds to a 
coverage of 84% of the solid angle in the center-of-mass 
frame for electrons, and 92% for muons. 

To increase the efficiency of the event selection, elec- 
tron candidate tracks are combined with photon can- 
didates to recover some of the energy lost through 
bremsstrahlung. In addition to the photon selection cri- 
teria listed above, photons used in bremsstrahlung re- 
covery are required to have A42 < 0.25. They are also 
required to be within 35 mrad in 9 from the track, and 
to have azimuthal angle intermediate between the ini- 
tial track direction and the centroid of the EMC cluster 
arising from the track. The initial track direction is es- 
timated by subtracting 50 mrad opposite to the bend 
direction from the (j) of the fitted track measured at the 
origin. The procedure increases the efficiency for recon- 
structing charmonium decays to e^e^ by about 30%. 



E. 



K? Candidate Selection 



We identify neutral hadrons through the presence of 
an energy deposition in the EMC or a cluster in the 
IFR. Neutral hadrons must be spatially separated from 
all tracks in the event. In reconstructing the decay 
B" -^ J/4>K1 neutral hadrons are taken as K^ can- 
didates, with requirements specifically tailored for this 
mode. 

Only the measured direction of the neutral hadron is 
used for K^ reconstruction, as its energy is poorly mea- 
sured. The direction of the K^ candidate is defined by 
the line joining the vertex of the J/ip candidate and the 
centroid of the EMC or IFR cluster. 

For a K'^ to reach the IFR it must traverse the EMC 
material, which amounts to approximately one nuclear 
interaction length. As a consequence, half of the K^ 



mesons undergo detectable interactions in the EMC. We 
consider EMC clusters with energy in the 0.2 - 2.0 GeV 
range. Most clusters arising from K^ interactions have 
energy below the upper bound; below the lower bound 
the contamination from noise becomes significant. All 
such EMC clusters which are spatially separated from a 
track are considered as K^ candidates, except those that 
combined with another neutral cluster give an invariant 
mass compatible with a 7r°. 

About 60% of K^ mesons from B'^ -^ J/ip K^ leave a 
detectable signal in the IFR. We select K^ candidates 
in the IFR starting with clusters of hits not spatially 
matched to a track. IFR clusters with hits only in the 
outer layers of the forward endcap are rejected to reduce 
the contribution from beam backgrounds. 



VI. EVENT SELECTION AND B MESON 
COUNTING 

A determination of B meson branching fractions de- 
pends upon an accurate measurement of the number of 
B mesons in the data sample. We find the number of BB 
pairs by comparing the rate of multihadron events in data 
taken on the T(4S') resonance to that in data taken off- 
resonance. The BB purity of the sample is enhanced by 
requiring the events to pass the following selection crite- 
ria, in which all tracks (including those that do not satisfy 
our usual selection requirements) in the fiducial region 
0.410 < 6 < 2.54 rad and all neutral clusters with energy 
greater than 30 MeV in the region 0.410 < 6 < 2.409 rad 
are considered: 

• The event must satisfy either the L3 DCH or L3 
EMC trigger. 

• There must be at least three tracks that satisfy the 
standard selection requirements in the fiducial re- 
gion. 

• The ratio of the second to the zeroth Fox- Wolfram 
moment urn must be less than 0.5. 

• The event vertex is calculated by an iterative pro- 
cedure that begins by considering every track in 
the event, and then discards those that contribute 
a large x^ to the fit (these are presumed to arise 
from the decay of long-lived particles) until the ver- 
tex fit is stable. This vertex must be within 0.5 cm 
of the beam spot center in xy and within 6 cm in 
z. The beam spot has an RMS width of about 120 
/im in X, 5.9 /im in y, and 0.9 cm in z. The point 
of closest approach of a high-momentum track to 
the beam spot is measured with a resolution of 23 
/im in X and y, and 29 /im in z, as determined with 
dimuon events. 

• The total energy of charged and neutral particles 
is required to be greater than 4.5 GeV. 



TABLE IV: Summary of observed invariant mass or mass dif- 
ference Am widths for all intermediate mesons considered in 
this paper. For most mesons the width is dominated by ex- 
perimental resolution, and the value reported in the table is 
the width a from a Gaussian fit to the data. For the K* 
modes the natural width of the resonance dominates, and the 
value reported is the full width of a Breit-Wigner fit to the 
data. The width for J/ip and ip{2S) decaying to e^e~ is 
greater than that for fi^ ij.~ due to the energy lost through 
bremsstrahlung. 



Quantity 


Decay mode 


Width (MeV/c^) 


J/iIj mass 


e+e- 


17 ±2 




/i+/i- 


13 ±1 


ip{2S) mass 


e+e- 


29 ±6 




M^M~ 


21 ±3 


Am(V'(25') - J/%Ij) ip{2S) -^ J/i> tt+tt" 


; 7±i 




j/ip -> i+r 




Am(xci- J/^) 


31^ -^ e+r 


14 ±1 


Ks mass 


^+^- 


3.5 ±0.2 




^° V' 


15 ±2 


K**^ mass 


K+ -K- and 


60 ±7 


K*^ mass 


K% 7r+ and 


50 ±10 



These requirements are 95.4 ± 1.4% efficient for BB 
events, as estimated from a Monte Carlo simulation. All 
events used in the branching fraction analyses are re- 
quired to pass this selection. 



VII. MESON CANDIDATE SELECTION 

The next step in the analysis is to combine sets of 
tracks and/or neutral clusters to form candidates for the 
initial or intermediate mesons in the decay. Our general 
strategy when forming these candidates is to assign the 
expected masses to tracks and neutral clusters, and to 
apply a vertex constraint before computing the invariant 
mass. In rare instances (less than 1% of all meson candi- 
dates) the vertex fit does not converge. The sum of the 
track and/or cluster four-vectors is used to compute the 
invariant mass for such candidates. If one or more decay 
products from a given particle are themselves interme- 
diate states, we constrain them to their known masses. 
At each step in the decay chain, we require that mesons 
have masses consistent with their assumed particle type. 
The mass resolutions observed for all of the intermediate 



mesons considered in this paper are listed in Table [V 



We choose meson selection criteria to maximize the ex- 
pected precision of our branching fraction measurements. 
Therefore we use well-understood quantities in our selec- 
tion, which lead to a smaller systematic uncertainty. We 
set the selection values to maximize the ratio S/^JS + B 
where S and B are the expected number of signal and 
background events respectively, as estimated from Monte 
Carlo. If a given mode has been previously observed, S 




ilMf mass (GeV/c ) 





- 


\ 






[L 


- 


100 


- 




- 




- (b) 


/ 


- 


5 






- 


^50 

in 
a 
W 






1 


25 


■^ 






- 




- J/ 


U : 


A 


[i.jTJi„,rAri^^'i-Wlr 


,S — „ - 



2.9 



3.1 



3.2 



J/v|/ mass (GeV/c ) 



FIG. 2: Invariant mass distribution for (a) J/ijj -^ e^e~ and 
(b) J/tp — > M^/i~ candidates mass observed in B*^ -^ Jlil)K% 
and B^ -^ J/tp K^ candidates passing the exclusive branch- 
ing fraction selection. The mass interval used to select J/'p 
candidates for B reconstruction is indicated by the arrows. 



is estimated using the known branching fraction. Oth- 
erwise, selection values similar to those in previously- 
observed modes are taken as a starting point, and then 
modified to reduce background (as measured in the kine- 
matic sidebands) or increase signal efficiency (as mea- 
sured using Monte Carlo simulated signal events). In 
most cases, we find that S/y/S + B does not change sig- 
nificantly when selection values are varied near their op- 
tima. This allows us to choose standard selection values 
across most final states. 



Charmonium Meson Candidate Selection 
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i J/tp Selection 



J/ip candidates are required to have an invariant mass 
in tlie range 2.95 < Mj/^ < 3.14GeV/c2 and 3.06 < 
Mj/^ < 3.14GeV/c2 for J/ip -^ e+e" and J/ip -^ ^J.+^l- 
decays respectively. Unless otherwise stated, for J/'ip -^ 
e^ e~ decays, one track is required to pass the tight elec- 
tron selection and the other the loose selection. Tracks 
not associated to an EMC cluster that pass the DCH-only 
selection are also accepted. For Jj-i^ -^ M"*"/^" decays, we 
require one track to pass the loose selection and the other 
to pass the MIP selection. 

The mass distribution for J/'0 candidates in the data 
is shown in Fig. 0. 



2 'tli{2S) Selection 



FIG. 3: Background-subtracted ^{2S) candidate mass and 
mass difference distributions observed in B'^ -^ ip{2S)Ks and 
B^ -^ il){2S)K'^ candidates passing the exclusive branching 
fraction selection, for (a) i/'(2S) -^ e'^e~ , (b) '0(25') -^ f^^ jj.~ , 
and (c) the tp{2S)-J/ip mass difference distribution for ^(25) 
—> J/4> 7r"'"7r~ The intervals used to select ^{2S) candidates 
for B reconstruction are indicated by the arrows. 
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'ip(2S) ^ /.t+/i candidates are required to have a 
mass within 50 MeV/c^ of the known ip{2S) value of 
3.69 GeV/c2 |§. For ^(25") ^ e+e" candidates the 
lower bound is relaxed to 250 MeV/c^ below the known 
value. For decays of the '0(25') to J/tp tt+tt^, the differ- 
ence in mass between the ip{2S) and J/ip candidates is 
required to be within 15 MeV/c^ of the expected value, 
and the tt+tt^ invariant mass m^+T^- is required to be be- 
tween 0.4 and 0.6 GeV/c^. The latter requirement takes 
advantage of the fact that m^+^- is most often in the 
upper portion of the kinematically allowed range P0[ . 
All ipi'2S) candidates are required to have a momentum 
in the center-of-mass frame between 1.0 and 1.6 GeV/c, 
consistent with B — > '4){2S)K decays. 

We have used the same lepton identification require- 
ments as for the J/ip reconstruction. These are applied 
either to the leptons from '0(25') — > £+£" decays, or to 
the leptons from the J/'0 in ^"(25) -^ J/ip 'n'^'K~ decays. 

The mass and mass difference distributions for tp{2S) 
candidates in the data are shown in Fig. || For Figures y, 
0, and o a background subtraction is performed using the 
observed distribution of candidates in the AE sidebands 



(see Section VII C) 



3 Xci Selection 



FIG. 4: Background-subtracted XcX-Jlip candidate mass dif- 
ference distribution observed in S" — > Xci^s and B^ — > 
Xc\K^ candidates passing the exclusive branching fraction 
selection. The mass difference interval used to select Xci can- 
didates for B reconstruction is indicated by the arrows. 



In reconstructing Xc\ —^ J/i^li J/i^ ELud photon candi- 
dates are selected as described above. The muon identi- 
fication requirements are subsequently tightened by de- 
manding that one lepton from the J/'0 pass the loose 
selection and the other the very loose selection (rather 
than the MIP selection). 

In addition, the photon cluster is required to satisfy 
E > 150 MeV and ^42 < 0.15 and to have a centroid 
in the angular range 0.41 < 9 < 2.409, excluding the 
forward direction due to the increased material (from 
electronics, cables, and final-focusing magnets) in that 
region. 
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FIG. 5: Ks candidate mass distribution observed in B" —* 
J/ip Kg candidates passing the exclusive branching fraction 
selection, for (a) K"^ -^ tt+tt" and (b) A's -^ vrV. The 
mass intervals used to select K^ — > 'K'^ir~ candidates for B 
reconstruction is indicated by the arrows in (a) ; the full range 
of (b) is used in selecting K^ -^ tt^'tt'' candidates. 



We require the mass difference between the recon- 
structed Xci and J/ip candidates to satisfy 0.35 < 
M^j/^ - Mj,^^ < 0.45 GeV/c2. 

The mass difference distribution for Xci candidates in 
the data is shown in Fig. 0. 



B. Light Meson Candidate Selection 

1 tt" -^ 77 Selection 

We reconstruct n^ candidates as pairs of photons. Indi- 
vidual photons separated by distances of 10 cm or more in 



FIG. 6: Background-subtracted (a) K*'^ and (b) K*'^ can- 
didate mass distributions observed in B" -^ J/tp K*" and 
B^ -^ J/ipK*'^ candidates passing the exclusive branching 
fraction selection. The mass interval used to select K* can- 
didates for B reconstruction is indicated by the arrows. 



the EMC are reconstructed as distinct clusters. Photons 
from 7r°'s with energies above 2 GeV can have less sepa- 
ration, in which case the two photons are reconstructed 
as a single cluster. We refer to these as "merged" tt^'s. 
They are distinguished from single photons based on their 
shower shape. 



Kl 



Selection 



We construct K^^ candidates from all pairs of oppo- 
sitely charged tracks, and retain those that have invari- 
ant mass between 489 and 507 MeV/c^ after applying a 
vertex constraint. To further reject background we ex- 
ploit the flight length of the K^ by demanding that the 
K^ vertex be more than 1 mm (in three dimensions) 
from the J/i/', tp{2S)^ or xa vertex. 

The mass distribution for K^ -^ tt+tt" candidates in 
the data is shown in Fig. ^. 



3 K° 



The K\ 







T°7r° -^ 47 decay chain is reconstructed 
from photon combinations satisfying E^ > 30MeV, 
E^o > 200 MeV and E^o > 800 MeV, with 110 < m^o < 
155MeV/c2 and 300 < ttiko < 800MeV/c2. We per- 
form a mass-constrained fit to each photon pair with the 
known ir'^ mass. This fit is repeated assuming different 
decay points along the iiTg flight path, as defined by the 
J/^l' vertex and the initial K*^ momentum vector direc- 
tion. The point where the product of the fit x^ proba- 
bilities for the two tt^'s is maximal is defined as the K'^ 
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decay vertex. K'^ candidates with flight length in the 
range from —10 to +40 cm are retained. 

We consider merged 7r° candidates with energy above 
1 GeV. If an EMC cluster candidate is identified as a 
merged tt'^ but can also be paired with another photon 
to form a 7r° candidate, we use the latter interpretation. 
Merged 7r°'s represent less than 10% of all tt^'s used in 
this analysis. 

The invariant mass of the K^ candidate at the opti- 
mal vertex point is required to lie in the range 470 to 
550MeV/c2. 

The mass distribution for K'^ -^ tt^tt*^ candidates in 
the data is shown in Fig. pi. 



4 K*'^ and K*'^ Reconstruction 

We reconstruct the K*'^ through its decays to K'^7r~ 
and K'^tt^ and the K*'^ through its decays to K'^tt'^ and 
A'+tt", where the K*^ is reconstructed in the tt+tt^ mode. 

tt'^'s are reconstructed from isolated photons and re- 
quired to have an invariant mass between 106 and 
153MeV/c^. If there is a K'^ in the final state we re- 
quire that the angle in the xy plane between the K^ 
momentum vector and the line joining the J/ip and K'^ 
vertices be less than 200mrad and that the K*^ vertex 
fit converge. 

In addition, for channels containing a tt" in the final 
state, we demand that the cosine of the angle 6k, mea- 
sured in the K* rest frame, between the kaon momentum 
and the K* direction as measured in the B frame be less 
than 0.95. 

All candidate K*'s are required to be within 
100MeV/c2 of the known K*° or K*+ mass ^. 

The mass distribution for K* candidates in the data is 
shown in Fig. ||. 

C. B Meson Candidate Selection 

B mesons are reconstructed by combining charmonium 
meson candidates with light meson candidates. Both 
the charmonium and light meson candidates are con- 
strained to their known masses, with the exception of 
K* candidates, for which the natural width dominates 
the experimental resolution. Two kinematic variables are 
used to isolate the B meson signal for all modes except 
S" — > J/ipKl. One is the difference between the recon- 
structed energy of the B candidate and the beam energy 
in the center-of-mass frame AE. The other is the beam 
energy substituted mass toes, defined as: 




FIG. 7: Helicity angles for the decay 7(45) ^ BB ^ J/^l) 
(e+e- or ii+ ^J.-)+Kl. 




TOES 



beam 



Pb 



(3) 



FIG. 8: Distributions of cos^^ observed in B'^ — > J/ip Ks and 
B^ ^ J/ip K^ candidates. The dashed histogram shows can- 
didates in the AE sideband. The solid histogram shows the 
distribution in the A£'-r7iES signal region, after subtracting 
the distribution observed in the sideband scaled by the ra- 
tio of signal to sideband areas. The normalization of both 
histograms has been set to unity. 



where pg is the momentum of the reconstructed B and 

frame 

of the run are taken into account when calculating tties 



beam i^ ^^^ beam energy, both in the center-of-mass 
The small variations of -E^cam over the duration 
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TABLE V: Definition of the signal region in \AE\ and mES for each mode used in this analysis. The rriES signal region is given 
in terms of |mES — "is|, where ms is 5279 MeV/c^. 



B decay 


Light meson 


Charmonium meson 


|A£;| (MeV) |mES-ms| (MeV/c^) 


mode 


decay mode 


decay mode 






B" -^J/i)K"s 


TT+n- 


e+e- 


34.5 


8.1 






/i+^^ 


29.0 


7.2 




7r%0 


e+e- 


100.0 


8.0 






m^m" 


100.0 


10.0 


B°~*J/t^Kl 


- 


e^e~ & l^'^f^~ 


10.0 


- 


B+ ^ J/iPK+ 


- 


e+e- 


38.4 


7.5 






/i+^- 


30.3 


6.9 


B" ^ J/iP K"" 


K+ Tv- 


e+e- 


30.9 


9.3 






/i+^- 


23.7 


8.1 




K°s7V° 


e+e- 


48.6 


12.0 






/i+^- 


45.6 


11.4 


B+ -^ J/i?K"+ 


Kl^- 


e+e- 


62.7 


7.2 






fj,+fi- 


20.4 


9.9 




T^+ttO 


e+e- 


85.2 


11.4 






fi+fi- 


50.1 


10.2 


B° ^ J/^jTv" 


77 


e+e" & ^+/i- 


112.0 


9.0 


B° -^ i>{2S)K°s 


TT 7r~ 


e e & e e vr tt 


28.0 


9.0 






f^^ IJ.~ & fJ-^ l^~ 7r+7r~ 


26.0 


9.0 


B+ -^ xp{2S)K+ 


TT 7r~ 
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Signal events will have AE close to and toes close to 
the B meson mass, 5.279 GeV/c^. 

We limit all our two dimensional plots in these vari- 
ables to the "signal neighborhood", defined by |Ai?|< 
Ai^inax and 5.2 < toes < 5-3 GeV/c^. For most chan- 
nels, A^max is 120 MeV, but for the B° -^ J/^PK°{K° 
— + TT*^ 7r°) and B'^ -^ J/'ipir^ channels, which have larger 
AE resolution, it is increased to 150 and 400 MeV respec- 
tively. We define the signal region by fitting the observed 
distribution of events in the signal neighborhood in toes 
and AE separately. In the fit, the signal component is 
modelled by a Gaussian, and the background component 
is modelled by an empirical phase-space distribution |2l[| 
(henceforth referred to as the ARGUS distribution) when 
fitting the ttjes distribution, or a polynomial when fitting 
the AE distribution. The ARGUS distribution is 



A(toes ; JTio , c) oc Toes v^l - (toes/^o)^) x 

exp(c(l - (toes/too)^)), (4) 

where toq is set to a typical beam energy and c is a fitted 
parameter. 

The widths of the fitted Gaussians provide a measure- 
ment of the resolution in AE and to-es, and the signal 
region is defined as ±3cr about the nominal value in each 
variable. The resolution in toes is typically 3 MeV/c^, 



and that in AE is typically 10 MeV for channels with no 
neutral particles in the final state and 30 MeV otherwise. 
The signal region for each mode is given in Table M. 

A somewhat different procedure is required for recon- 
structing B^ — > J/iJjKI, since the K^ energy is not mea- 
sured. Either the B mass or energy must be constrained, 
leaving only one independent variable. We choose to 
fix the B mass to its known value |1^ and plot the 
signal in the quantity A^^o= E^^^ + E*^,^ - E^^^^, 

where Ej,, is the energy of the mass-constrained J/^p , 
and £'^0 is the energy of the K^ as determined using 
the B mass constraint, both in the center-of-mass frame. 
AEj^o is a measure of the same quantity as AE; we 
use the different notation to reflect the fact that the B 
mass constraint is used only in this channel. For signal, 
AE^o is expected to peak at zero with a resolution of 
approximately 3.5 MeV. The signal region is defined as 
|A£;;^o| < 10 MeV. 



1 Helicity and Thrust Angle Definitions 

We use the helicity angles Ob and 6^ to help distinguish 
signal from background. 9b is the angle in the center-of- 
mass frame between the electron beam and B candidate 



14 



directions, and Oi is the angle in the charmonium meson 
rest frame between the i~ and hght meson candidate 
directions. Figure gives a schematic representation of 
these angles for the decay B'^ -^ J/ipK'^. 

The angle 9b has a sin^ 9b distribution for T(45) me- 
son decays. If X is a pseudoscalar {K'^, K^, 7r°) then 
the charmonium meson must be longitudinally polarized, 
and the resulting 9g distribution is proportional to sin 9^. 
If X is a vector {K* ) the decay angular distribution de- 
pends on more than one helicity amplitude. In this case 
the lepton angular distributions are not known a priori 
and must be experimentally determined. 

The B candidates formed from light quark back- 
grounds will generally follow a 1 -f cos^ 9b angular dis- 
tribution. The 9e helicity angle is especially useful in 
rejecting background since the distribution of cos 9i is 
peaked at ±1 for background and at zero for signal for 
modes where X is a pseudoscalar. As an example, the 
distribution of cos 9i observed in data for i?° -^ J/ip K^ 
and B^ — > J/ipK~^ candidates is shown in Fig. ^. 

For modes where the charmonium meson decays to 
more than two bodies, and 9i is therefore undefined, we 
suppress backgrounds using the thrust angle 9t, defined 
as the angle between the thrust axis of the reconstructed 
B and that of the rest of the event in the center-of-mass 
frame. We use the conventional definition of the thrust 
axis for a collection of particles as the direction about 
which the transverse momenta of the particles is min- 
imized. In BB events cos St is uniformly distributed, 
whereas in continuum background events 9x tends to 
peak at tt radians due to the two-jet nature of these 
events. Hence 9t can be used to discriminate against 
background in modes where the helicity angle is not ap- 
plicable. 

The helicity and thrust angle values used to select can- 
didates are listed in the appropriate exclusive reconstruc- 
tion and selection subsections in this paper. 



2 Multiple Candidates 

We only allow one exclusive candidate per event in a 
given decay mode. In the cases where we have multi- 
ple candidates (less than 10% of all events with a candi- 
date for most modes, but up to 30% for the K* modes 
which have significant crossfeed among decay channels), 
the candidate with the lowest |A£'| is taken over all oth- 
ers. The only exception is in the i?° -^ J/ipK^ selection, 
where we choose the candidate with the largest K^ en- 
ergy as measured by the EMC. If none of the candidate 
K^ mesons have EMC information, we choose the candi- 
date that has the largest number of layers with hits in the 
IFR. These criteria are used because background candi- 
dates tend to have low EMC energy or IFR multiplicity. 
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FIG. 9: Signals for B° ~ 
K"s -^ tt" tt") and (c) B" 
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AE (MeV) 



J/ipK^ ((a) A'g -^ TT+TT" and (b) 
-» J/ip K1 . In (a) and (b) the upper 
plots show the distribution of events in the AE-rriES plane, 
and the lower plots show the distribution in tties of events in 
the signal region in AE. In (c) the points are the data, the 
dot-dashed line shows the Monte Carlo simulated distribution 
of background events which include a real J/ip , the dashed line 
shows the model for the total background, where the non- J/?/) 
component is taken from the J/tp sidebands in data, and the 
solid line shows the sum of the background and signal Monte 
Carlo models. 
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FIG. 10: Signal for (a) B+ -^ J/ipK+, (b) B° -^ J/V'tt", (c-d) B° -^ J/ipK*°, and (e-f) B+ ^ J/'ipK*+. The upper plots 
show the distribution of events in the AS-jties plane, and the lower plots show the distribution in mES of events in the signal 
region in A_E. 



3 B° ~* J/tpKl:(n+n-) 

All combinations of J/tp and K'^ -^ tt+tt^ candidates 
are used to form B candidates. We require the absolute 
value of cos 9i to be less than 0.8 and 0.9 for J/ip -^ e^e~ 
and J/'0 — > M^/^^ events respectively. The distribution 
of the selected candidates in Ai? and toes is shown in 
Fig. I (a). 



A B° 



J/^7^°r^°7rO; 



All combinations of J/t/i and Xg -^ tt^tt" candidates 
are considered. For J/i/; — + e^e^ candidates, one track is 
required to pass the tight or DCH-only selection, and no 
particle identification requirement is placed on the sec- 
ond track. The mass-constrained Jji]} vertex is assumed 
to be the production point of the i^g. We require that 
the absolute value of cos0£ be less than 0.7 and 0.8 for 



Jji]} — *■ e^ e" and J/f/' — > [i-^ [i-^ events respectively. The 
distribution of the selected candidates in A£' and toes is 
shown in Fig. || (b). 



B" 



JI^Kl 



Since most of the background in this mode arises from 
-B decays that include charmonium mesons, we reject 
events if they contain a candidate for B^ — > J/ipK^;: 
B+ -^ J/tPK+, B° -^ J/^PK*^, or B+ -^ J/^IjK*+. The 
decay modes used to resconstruct these candidates are 
the same as those used in the branching fraction analysis 
for each mode, but the selection criteria are loosened. 

Within the remaining events, we select J/tp candidates 
using a procedure that differs slightly from the standard 
selection. A vertex constraint is applied, and only can- 
didates for which the fit converges are retained. In ad- 
dition the momentum of the J/tp in the center of mass 
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FIG. 11: Signal for (a) B° ^ tp{2S)K'i and (b) B+ -^ 
'<p{2S)K^ . The upper plots show the distribution of events in 
the A_B-mES plane, and the lower plots show the distribution 
in niES of events in the signal region in AE. 



frame is required to be between 1.4 and 2.0 GeV/c, con- 
sistent with i?'^ — > J/^K^ decays. In the e'^e~ mode, 
one electron candidate is required to pass the very tight 
selection and the other the loose selection, and the J/ip 
mass is required to be between 3.00 and 3.13 GeV/c^. For 
the ^1^ n'' mode one muon candidate must pass the tight 
selection and the other the loose selection, and the J/i/j 
mass is required to be between 3.06 and 3.13 GeV/c^. 

We consider all pairs of K^ and J/tp candidates, as 
described above, as candidates for i? ^- -0 K^ decays. We 
then construct the quantity AEj^o described previously. 

For candidates containing a K^ that is identified in the 
EMC, we require that the transverse missing momentum 
be consistent with the momentum of the K^ candidate 
calculated from the B mass constraint. We compute the 
missing momentum from all tracks and EMC clusters, 
omitting the K^ candidate cluster. This quantity is then 
projected along the direction of the K^ candidate in the 
plane transverse to the beam. Studies of B'-^ -^ J/tp K^ 
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FIG. 12: Signal for (a) B° ^ Xci^s, (b) B+ ~* XciK+, and 
(c) B -^ X<:^K* . The upper plots show the distribution of 
events in the A_E-r7iES plane, and the lower plots show the 
distribution in niES of events in the signal region in A_B. 
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events in the simulation imply that the event missing mo- 
mentum should be equal to the calculated momentum of 
the if", with a resolution of 0.30 GeV/c. Therefore, we 
select events where the total missing momentum is not 
less than 0.65 GeV/c below the calculated K'^ momen- 
tum. The missing momentum requirement is not applied 
when the K^ candidate is identified in the IFR, since the 
background is much lower in this sample. 

For all events, we use the angles 9b and 0i to suppress 
background. We require that | cos^bI and | cos9e\ be less 
than 0.9. To further reduce background, we also demand 
that I cos 6'_B I + I cos 0^1 be less than 1.3. 

The distribution of the selected candidates in AEf^o 
is shown in Fig. g (c). 



6 B^ 



J/tl}K+ 



the ip{2S), cosOt is required to have an absolute value of 
less than 0.9. The distribution of the selected candidates 
in AE and ttt-es is shown in Fig. |ll|. 

10 B° -^ XciK°s and B+ ~* XciK+ 

B^ —> Xci-^'s candidates are formed by combining 
mass-constrained Xci candidates with mass-constrained 
K^l — > 7r+7r~ candidates. 

K^ candidates are defined as tracks which lie within 
the angular range 0.35 < 9 < 2.5 rad. These are 
combined with mass-constrained Xci candidates to form 
_B+ -^ XciK^ candidates. 

The cosine of the 9t is required to have absolute value 
less than 0.9. The distributions of the selected candidates 
in AE and mES are shown in Fig. O (a-b). 



Every combination of a J/^p candidate and a track is 
considered. We require | cos 9e\ to be less than 0.8 and 0.9 
for J/i/i —f e+e~ and J/ip — > fJ-'^fJ.~ events respectively. 
The distribution of the selected candidates in AE and 
Toes is shown in Fig. ^ (a). 



7 B" 



J/lpTT° 



For J/ip — > fi^fj.^ the standard selection is tightened 
by requiring that one charged track satisfy the very tight 
criteria and the other the loose criteria. Only 7r° 's formed 
from isolated photon pairs with mass between 120 and 
150 MeV/c^ are considered. 

The absolute value of cos 9t is required to be less than 
0.95. Since background events are correlated in 9t and 
9e, we also demand that |cos6't| + | cos 6*^1 be less than 
1.8. The distribution of the selected candidates in AE 
and Toes is shown in Fig. hO (b). 



8 B°^J/^K* 



B^ 



j/i>K*^ 



The S" is reconstructed from pairs of J/4> and if *° 
candidates, while the B^ uses J/ip and K'^ candidates. 
We further require that both J/-0 daughter leptons satisfy 
either the loose muon selection criteria or tight electron 
selection criteria. 

The distribution of the selected candidates in AE and 
Toes are shown in Fig. |l^ (c-f). 



9 B° ~* ^{2S)K°s and B+ -^ iP{2S)K+ 

Charged B candidates are formed from the combina- 
tion of a ■0(25') candidate with a track, and neutral can- 
didates from the combination of V'(2S') and K^ -^ tt+tt^ 
candidates. 

In the leptonic decay mode of the tp{2S), |cos6'f| is 
required to be less than 0.8. In the J/f/; decay mode of 



11 B°^XciK*° 

B candidates are reconstructed by combining mass- 
constrained Xci candidates with if *° candidates recon- 
structed in the if+7r~ mode. We require that the if+ 
candidate be inconsistent with a pion hypothesis, using 
the combined information from dE/dx measured in the 
SVT and DCH and Cherenkov angle measured in the 
DIRC. We require both tracks from the J/ip to pass ei- 
ther the tight electron selection or the loose muon se- 
lection. Xci candidates are selected if the mass differ- 
ence between the Xci and the J/ip lies between 0.37 and 
0.45 GeV/c^. if *" candidates are reconstructed using the 
standard procedure, and are accepted if the if *° mass is 
within 75 MeV/c^ of the known value ||l^ . 

The distribution of the selected candidates in AE and 



TO.ES is shown in Fig. 12 (c). 



VIII. BACKGROUND ESTIMATION 

Backgrounds to the decay modes we measure arise pre- 
dominantly from three sources: other B decays that in- 
clude charmonium mesons in the final state, B decays 
without charmonium mesons, and light quark events. 
Monte Carlo simulation studies verify that for B decays 
without charmonium mesons and for continuum events, 
B candidates follow the ARGUS distribution in toes. 
On the other hand, the background from inclusive char- 
monium decays includes modes that are kinematically 
very similar to the signal modes, which means that their 
distribution in toes may have a peak in the signal re- 
gion. As an example. Fig. O shows the distribution in 
AE and toes for signal and background events satisfying 
the B~^ -^ XciK^ selection requirements. It is critical 
that the so-called "peaking background" from other J/-0 
modes be well understood, since it contributes directly 
as a correction to the fitted number of signal events in 
the signal band. 



18 



80 
60 


1 ' 1 ' 
Signal + background 
tai Other backgrounds 

1 1 J/\|/ background 

:(a) 


1 1 


20 






- 


(\ 


ife 


\%^^^f»r^ 


-+ : 

^ - 



Combinatorial 



Peaking 



5.2 5.225 5.25 5.275 5.3 
mj53 (GeV/c^) 




-0.1 -0.05 0.05 0.1 

AE (GeV) 



FIG. 13: Distribution in (a) tties and (b) AE of candidates 
for B'^ -^ XciK'^ . The points are the data, the shaded his- 
tograms are Monte Carlo simulated background events, bro- 
ken down into the combinatorial and inclusive J/ip contri- 
butions, and the open histograms are the sum of the Monte 
Carlo simulated signal and background distributions. The 
Monte Carlo distributions are normalized according to the 
equivalent luminosity of the samples. In (b) the AE signal 
region lies between the solid arrows, and the sideband region 
in which we compare the observed peaking background to 
the Monte Carlo prediction lies outside of the dashed arrows. 
Note that the inclusive J/tp background peaks in the signal re- 
gion of r7iES, but that neither background peaks in the signal 
region of A_B. 
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FIG. 14: Difference between the predicted and observed lev- 
els of background, divided by the combined statistical error 
from data and Monte Carlo simulation. The comparison of 
combinatorial backgrounds is done in the signal region, while 
for peaking backgrounds the AE sideband region is used. For 
the J/ipTT^ mode the value shown is the sum of the e^e~ and 
fi'^ IJ.~ modes. 

TABLE VI: Dominant sources of background in the decay 
modes we consider, along with the fraction of the total back- 
ground due to the dominant source. These fractions have 
substantial uncertainty due to the limited statistics of the 
available Monte Carlo simulation sample. 



Channel 



B" 



B+ 

B°- 

B°- 

B+ 

B°- 

B+ 

B"- 

B+ 

B«- 



J/ijjK" 



Dominant % of total 

background 






> J/ipK+ 

J/iPK*° 
*.J/xIjK*+ 

ip{2S)K° 
» iP{2S)K+ 

* XciK+ 

X.iK*° 



TT'*'Tr~ Charmonium 70 

tt" tt" Continuum qq 50 

Charmonium 90 

Charmonium 50 

Continuum qq 55 

Charmonium 90 

Charmonium 85 

Charmonium 60 

Charmonium 50 

Charmonium 95 

Charmonium 75 

Charmonium 90 
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For all modes except i?" -^ J/tpKl, we estimate the 
m.agnitude of the backgrounds by using Monte Carlo sim- 
ulation, off-resonance data, and mass sidebands for J/ip 
or 4'(2S) candidates in on-resonance data. The avail- 
able Monte Carlo samples are 10 million BB decays, the 
equivalent of 8 fb^ of continuum events, and the equiv- 
alent of several times our data sample of inclusive B to 
charmonium decays. 

We compare the predicted and observed levels of back- 
ground in two regions of the AE-m^s plane: the AE 
sideband, defined as that part of the signal neighbor- 
hood sufficiently far from the signal region in \AE\ that 
it contains a negligible amount of signal (typically 4cr 
from zero, though for modes with a tt" in the final state 
this is reduced to 3cr), and the signal region. 

In each region, we fit a Gaussian and an ARGUS back- 
ground distribution to the observed toes distribution of 
B candidates in data and Monte Carlo samples. In the 
AE sideband the integral of the Gaussian distribution 
across the to.es signal region provides an estimate of the 
peaking background. In the AE signal region the integral 
of the ARGUS background function across the toes signal 
region provides an estimate of the combinatorial back- 
ground. A comparison between data and Monte Carlo 
simulation of the fitted results for the combinatorial and 
peaking background components is displayed in Fig. |lj. 
In most cases, the predicted and observed backgrounds 
are in good agreement, within the statistical errors. Dis- 
crepancies in the predicted and observed levels of peaking 
backgrounds in the AE sideband region are accounted for 
in our estimation of systematic uncertainties. 

For the B'-* — > J/ipK^ sample, we estimate the mag- 
nitude of the background by performing a binned log- 
likelihood fit to the AE}^o distribution in the range —20 
to 80 MeV. This fit is described in detail in Section ^. 
The shapes of the signal and inclusive charmonium back- 
ground components are taken from Monte Carlo simu- 
lation. The shape of the non-charmonium background 
component is taken from an ARGUS fit to the AEj^o 
distribution for events in the J/tp mass sideband. To con- 
strain the magnitude of this last component, we first esti- 
mate the fraction of non- J/ip candidates in the J/tp mass 
window relative to the mass sideband for events with ar- 
bitrary AEj^o . We then scale the number of events with 
AEj^o between —20 and 80 MeV that also have a dilep- 
ton invariant mass in the J/'0 sideband region by this 
fraction to determine the expected number of candidates 
arising from non-charmonium backgrounds. 

The dominant source of background for each mode we 



consider is listed in Table VI 



IX. EFFICIENCY CALCULATION 

The selection efficiencies for each mode are obtained 
from detailed Monte Carlo simulations, in which the de- 
tector response is simulated using the GEANT3 |g2[ pro- 
gram. In addition, we have used the data where possible 



to determine the detector performance. 

We have determined the efficiency for identifying lep- 
tons with the sample of inclusively produced J/ip 's in the 
data. J/ip 's are selected by requiring that one track pass 
the very tight electron or muon selection, with no lepton 
identification requirement placed on the other track (the 
test track). The fraction of test tracks that satisfy a given 
lepton selection provides a measure of the efficiency for 
that selection. 

We have determined the track finding efficiency from 
multihadron events in the data. For the standard track 
selection, the fact that the SVT is an independent track- 
ing device allows precise determination of the DCH ef- 
ficiency by observing the fraction of tracks in the SVT 
that are also found in the DCH. For low-momentum pi- 
ous, such as those produced in the decay ■ip{2S) — > J/4> 
TT+TT^, D* decays are used to provide information about 
the efficiency as a function of momentum. This measure- 
ment takes advantage of the correlation between the pion 
helicity angle in the D* rest frame and its momentum in 
the center-of-mass frame. Since the helicity angle distri- 
bution is known, any deviation between the expected and 
observed distributions can be interpreted as arising from 
a momentum dependence in the track reconstruction effi- 
ciency. In addition, the efficiency for reconstructing a K'^ 
~f TT+TT^ decay has been determined as a function of the 
if g flight length from studies of inclusive ATg production 
in the data. 

The efficiency for detecting photon clusters has been 
determined from the data with a control sample of two- 
prong T'^T~ events. In the subsample of events tagged 
by a leptonic decay of one of the taus, we compare the 
number of events with one or two neutral pions, and one 
charged pion, from the second tau decay. The ratio of 
these two branching fractions is known to a precision of 
1.6% ||l^. By comparing data with simulation, we de- 
termine a correction factor to be applied to the photon 
identification efficiency. This factor is found to be inde- 
pendent of the photon energy. 

Both the J/4> mass distribution and AE signal distri- 
bution in the 5+ -^ J/tpK'^ sample have better resolu- 
tion in the simulation than in the data, indicating that 
the track px resolution in the simulation is overestimated. 
To account for this, we degrade the px resolution of the 
simulated tracks by amount chosen to bring the simulated 
J/ip mass and AE mass distributions into agreement with 
those observed in data. 

We measure the efficiency of the EMC and the IFR to 
detect a K^ candidate cluster using a control sample of 
e+e" -^ $7, $ -^ K°K° events. 

The efficiencies of the ir'^ veto and missing transverse 
momentum requirements applied for K^ reconstruction 
in the EMC were determined using 5+ -^ J/tp K~^ events. 

The AEj^o distribution for simulated events is ad- 
justed slightly to account for differences between data 
and Monte Carlo simulation in the beam energy spread 
and K'l angular resolution. The correction to the beam 
energy spread is derived from a study of 5+ -^ J/ipK^ 
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events, and the adjustment for the K^ angular resolution 
is determined with the e+e~ -^ $7 control sample. 

The combination of these effects requires a correction 
factor to be applied to the efficiency determined from the 
Monte Carlo simulation. The size of the correction varies 
among decay modes, and is at most 16%. 



X. BRANCHING FRACTION 
DETERMINATION 

To derive branching fractions we have used the sec- 
ondary branching fractions S published in Ref. |]19[. An 
exception to this is the branching fraction of ip{2S) — > 
£~^i~, where we have used our recent measurement of 
(6.6 ± 1.1) X 10-3 m for the ij{2S) -^ fi+fi' mode and 
the measurement from E835 p4] for the ^(25) — > e~^e~ 
mode. These measurements are more recent and more 
accurate than those included in Ref. jl9| . 

We have assumed that T{4S) decays produce an equal 
mixture of charged and neutral B mesons. The depen- 
dence of our results on this assumption is included in 
Section ^. 

For all modes except B° -^ J/ipK^, B" -^ J/ipK*" and 
_B+ -^ J/tpK*^, the number of signal events Ng within 
the signal region of the A_E-toes plane is determined from 
the observed number of events after background subtrac- 
tion. The background has two components, as described 
in Section VII] : a combinatorial component, which is 



obtained by integrating the fitted ARGUS distribution 
in the signal region, and a peaking component that is 
obtained from inclusive B -^ J/tpX simulation after re- 
moving the signal channel. The procedure is illustrated 
in Fig. |l|. 

We determine the branching fraction by dividing Ng 
by the selection efficiency e, S, and the number of BB 
events in the sample N^g. Where possible, the branch- 
ing fraction is determined independently for the different 
secondary decay modes, and the results combined statis- 
tically, taking into account correlations in the systematic 
errors. For the channels that are statistically limited, 
we determine the branching fraction using the combined 
sample of candidate B events, irrespective of the sec- 
ondary decay mode: 



B 



E.Ns^^ 



NBBT.^^^Si 



(5) 



60 



^40 



in 
I 20 



'>^*4^iMi^*^ 



5.2 5.225 5.25 5.275 
mg3 (GeV/c^) 



5.3 



FIG. 15: Distribution in mss of candidates for B"*" -^ 
XciK^ , with the ARGUS and Gaussian fit superimposed. 
The number of signal events is calculated by counting the 
events in the signal region of jties (marked by arrows) and 
subtracting the integral of the fit ARGUS function across this 
region (the shaded portion of the fit) and the peaking contri- 
bution from inclusive J/i}} backgrounds, as shown in Fig. |l3|. 



form: 



/;(-B°,6+) = Jl 



Mi, e 



N^J\ 



(6) 



where i represents a decay mode of the K* (to K^ or 
K+), j represents either the i?" — > J/'ipK*° or 5+ -^ 
J/ipK*~^ mode, N is the observed number of events in 
the signal region, and /i is the expected number of events. 
The last is given by: 



i'j' 



iji'j 



,S,,yN^ 



BB 



(7) 



where the sum is over all decay modes considered. 

The branching fractions for the S" -^ Jj^jK*^ {&) 
and B^ -^ J/ipK*~^ (13^) niodes are determined simul- 
taneously from a likelihood fit, which is required to ac- 
count for the cross-feed between the K* decay channels. 
The cross-feed is largest for the mode B^ -^ J/ipK*^ , 
where the K*+ decays to ivT+Tr". In this case, 12% of the mass sidebands of the J/t/j in the z*'' bin, and dj, the 



where N}, is the number of background events estimated 
in the same manner as for the other channels. The four 
indices attached to the selection efficiencies denote the 
fraction of events in the i' j' mode that pass the ij selec- 
tion requirements, as determined with the Monte Carlo 
simulation. 

We determine the number of signal and background 
events for the B'^ — > J/ij^K^ decay mode by performing 
a binned likelihood fit to the AE-l-o distribution. The fit 
takes as input a^, the fraction of simulated B^ -^ J/'4'K^ 
events in the i^^ bin, 6^, the fraction of simulated inclusive 



charmonium background events in the 



„th 



the 



fraction of non-charmonium background events from the 



selected candidates arise from other B^ J/ip K* decays. 
The likelihood function includes the cross-feed contribu- 
tions as well as all other background sources, and has the 



number of data events in the 



■th 



bin. The likelihood 
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TABLE VII: Breakdown of contributions to the systematic errors. Included are the contributions from the secondary branching 
fractions (S), lepton identification efficiency (PID), track pt resolution (Trk pr), track and K^ -^ -k^-k~ reconstruction 
efficiency {e{Trk+Ks)), photon identification efficiency (e(7)), background determination (BGR), Monte Carlo statistics (A^sim) 
and selection requirement variation (Sel. var.). The 1.6% error from the determination of the number of BB events, which is 
common to all modes, is not listed but is included in the total. In addition, the statistical uncertainty is shown. All values are 
expressed relative to the measured branching fraction, in percent. 



Channel 




5 


PID 


Trk Pt 


e{Trk+K'>) 


^(7) 


BGR Nsim 


Sel. var. 


Total Stat, error 


B" -^ J/i^K" 


K^s -^ ^+^- 


1.7 


1.3 


0.9 


5.5 


- 


1.1 


1.3 


3.5 


7.3 


6.4 




K^ ^ ^« ttO 


1.7 


0.5 


0.1 


2.4 


5.0 


2.0 


1.6 


2.5 


7.0 


15.2 


B+ ^ J/ijK+ 




1.7 


1.4 


1.0 


3.6 


- 


1.0 


0.8 


2.2 


5.3 


3.1 


B«-. J/V-tt" 




1.7 


2.5 


0.4 


2.4 


2.5 


1.7 


1.1 


10.0 


11.3 


32.7 


B" ^JIi>K"' 




1.7 


1.3 


0.8 


4.7 


0.2 


1.4 


0.2 


4.0 


6.9 


4.0 


B+ ~* J/ijK'+ 




1.7 


1.3 


1.1 


4.9 


1.2 


2.9 


0.1 


5.0 


8.2 


6.6 


B" -> i>{2S)K'' 




9.6 


1.0 


1.3 


7.9 


- 


4.8 


1.4 


8.5 


15.9 


15.4 


B+ -> i^{2S)K+ 




9.6 


1.0 


1.3 


5.8 


- 


1.3 


1.6 


3.7 


12.1 


8.1 


dO ^ ,^ rv-O 




6.2 


2.4 


1.2 


5.6 


1.3 


14.5 


2.2 


13.2 


22.0 


25.1 


B+ ^ XciK+ 




6.1 


2.6 


0.5 


3.6 


1.8 


3.8 


2.4 


5.3 


10.6 


10.0 


B -* Xci^ 




6.2 


2.4 


0.8 


4.8 


2.7 


14.3 


1.8 


8.1 


18.7 


28.8 



TABLE VIII: Breakdown of contributions to the systematic 
error for the B —^ J/ipKi^ analysis. The statistical error is 
also shown, with all values expressed relative to the measured 
branching fraction, in percent. 



Source 



Uncertainty 



Tracking efficiency 2.4 

Lepton identification efficiency 1.2 

J/i/) mass requirement efficiency 1.3 

Kl efficiency 9 

n" veto efficiency 0.7 

Missing momentum requirement efficiency 0.5 

Beam energy scale (spread) 1.0 (3.0) 

Kl angular resolution 4 

Branching fractions for B^ J/V> X 3.8 

non-J/i/) background shape 2 

Simulation statistics 2.2 

Secondary branching fractions 1.2 

Number of BB events 1.6 



Total 



12.0 



Statistical error 



12.0 
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function has the form: 



/:(7v„7v^x,iv„o„-v>) 



JVbin 

n 


r/,1 


Mi 






1=1 












-(« 


aon- 


^ 


-M)2 


X — 


eK'^ 


+ iV 


„o„ 


-.) 



(8) 



^27r (ct2 + iVnon-^) 

where -/V^x is the number of inclusive charmonium back- 
ground events, A'non-^ is the number of non-charmonium 
background events, M is the expected number of non- 
charmonium background events determined from the 
mass sidebands of the J/ip , a is the uncertainty on Af, 
and fii is the expected number of events in the i'^ bin, 
defined as: 



Hi = NsQi + N^xbi + Nnon-i,Ci 



XI. SYSTEMATIC ERRORS 



(9) 



Systematic errors on the results arise from the un- 
certainty on the number of BB events, the secondary 
branching fractions of the modes considered, the esti- 
mate of the selection efficiency and the knowledge of the 
background level. The size of the various contributions 
to the systematic error, expressed as a fraction of the 
branching fraction value , is listed for all modes except 
B° -^ J/^pK° in Tabic ^ and for the B° -^ J/ipKl 



mode in Table VIII 



The uncertainty on the number of BB events intro- 
duces a systematic error of 1.6% in common for all modes. 
The uncertainties in the branching fractions of the sec- 
ondary decay modes lead to a systematic error of between 
1.7% and 9.8%, depending on the mode considered. 

The systematic error due to the finite size of the avail- 
able Monte Carlo sample is between 0.1% and 2.4% for 
the different modes. 

We have determined the efficiency for a charged parti- 
cle to be reconstructed as a track that passes the stan- 
dard track selection to a precision of 1.2% per track. The 
uncertainty in the reconstruction efficiency for the low- 
momentum pions from the ipC^S) — > J/xJj7r^n^ decay is 
determined to be 2% per track. The systematic error 
associated with reconstructing a K'^ -^ tt+tt" decay has 
two sources: knowledge of the reconstruction efficiency 
for the two tt tracks, and differences in the selection cri- 
teria efficiencies observed between the inclusive K'^ data 
and the Monte Carlo simulation. The observed discrep- 
ancies and their statistical uncertainties are summed in 
quadrature to yield a systematic error of approximately 
5%. 

The systematic error on lepton identification efficien- 
cies arise from the statistics of the inclusive J/ip sam- 
ple, and from comparing the efficiencies in different low- 
multiplicity control samples. It varies from 0.5% to 2.8% 
per J/ip or ip{2S) depending on the criteria used to select 
the leptons. 



The quality of the simulation of photon detection and 
energy measurement in the EMC has been validated by 
a detailed comparison between real and simulated data. 
In particular, the position and resolution of the tt^ and ij 
mass peaks in the photon pair mass spectrum has been 
compared as function of photon energy, calorimeter oc- 
cupancy and time of data collection. The agreement in 
terms of energy scale is found to be better than 0.75% in 
all cases; energy resolution is also well described at the 
level of 1.5%. The absolute photon detection efSciency 
is known to 1.25%. The resulting systematic errors on 
the branching fractions are in the range of 1.3% to 5% 
depending on the decay mode. 

We account for the uncertainty in the pt resolution 
by varying the amount by which the Monte Carlo simu- 
lated momentum resolution is degraded within the range 
in which the data and Monte Carlo J/tp mass and AE 
widths are compatible. The observed variation in selec- 
tion efficiency is between 0.1% and 1.3%. To account 
for the possibility that other variables used in selecting 
candidates may not be perfectly modelled in the simula- 
tion, we vary the selection requirements and repeat the 
branching fraction measurement. In most cases the range 
of variation is ±lcr, where a is the width observed in data 
for the variable under consideration, while for helicity 
angles a variation of ±0.05 in their cosine is used. The 
observed variations in the results are between 2.5% and 
14.1%. Modes with a K* in the final state merit special 
mention, since there can be some variation of selection 
efficiency with the polarization of the vector meson, and 
the polarization amplitudes are subject to experimental 
uncertainty. The Monte Carlo simulation from which we 
derive our efficiency assumes the polarization amplitudes 
measured by CLEO |2^. We have studied the changes 
in efficiency that occur when the amplitudes are var- 
ied by twice the difference between the values measured 
at CLEO and BABAR ^. We find that these changes 
are consistent with those observed when the selection re- 
quirement on 6k is varied. 

For the B^ -^ J/ipKl analysis, we include additional 
systematic errors associated with the selection efficiency. 
These originate from the uncertainty in the K^ recon- 
struction efficiency and angular resolution determined 
from data, the knowledge of the absolute scale and spread 
of the beam energy, and from the various selection re- 
quirements used to isolate the signal. 

Another systematic error arises from our knowledge of 
the backgrounds. For all modes except S" -^ J/ipK^, 
we use data in the AE sideband to estimate this un- 
certainty. We determine the uncertainty in the size of 
the combinatorial background by repeating the fit to the 
data with the shape of the ARGUS function (the param- 
eter c in Equation 0) fixed to the value obtained from 
fitting the AE sidebands, allowing only the normaliza- 
tion to vary. This accounts for any correlation between 
the ARGUS and Gaussian fits in the AE signal region. 
We estimate the uncertainty in the predicted size of the 
peaking background by comparing the observed Gaussian 
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TABLE IX: Measured branching fractions for exclusive de- 
cays of B mesons involving charmonium. The first error is 
statistical and the second systematic. 



Channel 




Branching fraction/10 * 


B" -* J/^K" 


K]i -. n^Tv- 


8.5 


± 


0.5 


± 


0.6 




Kl ^ ^« TfO 


9.6 


± 


1.5 


± 


0.7 




Kl 


6.8 


± 


0.8 


± 


0.8 




All 


8.3 


± 


0.4 


± 


0.5 


B+ -* J/ijK+ 




10.1 


± 


0.3 


± 


0.5 


fiO^ J/^TfO 




0.20 


± 


0.06 


± 


0.02 


B° ^ J/i^K*° 




12.4 


± 


0.5 


± 


0.9 


B+ -^ J/ipK*+ 




13.7 


± 


0.9 


± 


1.1 


B° ~* i}}{2S)K° 




6.9 


± 


1.1 


± 


1.1 


B+ -» i>{2S)K+ 




6.4 


± 


0.5 


± 


0.8 


B» - XciK'> 




5.4 


± 


1.4 


± 


1.1 


B+ ^ XciK+ 




7.5 


± 


0.8 


± 


0.8 


B" - XciK*'> 




4.8 


± 


1.4 


± 


0.9 



component in the A_E sideband to that estimated from 
the inclusive B -^ J/tjjX simulation. This procedure 
takes advantage of the fact that the distribution of can- 
didates from this background in AiJ depends primarily 
on kinematics rather than the poorly-known composition 
of the background. In particular, the background does 
not peak in the signal region of /S.E (see Fig. |l^, which 
implies that the relative normalization observed in the 
/S.E sideband can also be expected to hold in the signal 
region. The systematic error attributed to the knowl- 
edge of the backgrounds varies from 1.0% to 14.5% for 
the various modes. In addition, for the B^ -^ J/ipK*'^, 
5+ -^ J/ifjK*^ and B^ -^ XciK*'^ modes, a system- 
atic error is included to account for the uncertainty in 
the non-resonant B — > J/iliKir branching fractions, and 
the contribution of feed-down from higher K* resonances. 
This ranges from 1.4% to 3.7% depending on the mode. 
For the B^ -^ J/ipK^ decay mode, we determine the 
uncertainty arising from knowledge of the shape of the 
non- J/iJj background both by changing the fitted param- 
eters of the ARGUS function for this background compo- 
nent by one standard deviation and also directly in the fit 
by using the AEp^o distribution from the non-J/'ijj events 
in the data. The analysis is also repeated after varying 
the values of the branching fractions for the component 
modes in the simulation of B ~* J/ip X decays by the un- 
certainty quoted in Ref . [|l9| . This is done separately for 
the main background modes and then for all the remain- 
ing modes together. Since the non-resonant B — > J/ipK-K 
component is poorly measured, we vary it in the range 
from -50% to -f 400%. 



XII. RESULTS 



In Table IX we summarize our branching fraction mea- 
surements. The observed number of events in the signal 
region, the predicted background, and the selection effi- 



ciency are given in Table ^. 

From these results, we have determined the following 
ratios of charged to neutral branching fractions, where 
the first error is statistical and the second systematic: 



B{B+ ^ J/ijK+) 

B{B° ^ J/i^K^) 

B{B+ ^ J/jj K*+) 

B{B° ^ J/ipK*o) 

B{B+ -^ iP{2S)K+) 

B{BO -^ i;{2S)K0) 

B{B+ ^ XciK+) 



= 1.20 ±0.07 ±0.04 (10) 

= 1.10 ±0.09 ±0.08 (11) 

= 0.94 ±0.16 ±0.10 (12) 

= 1.39 ±0.37 ±0.22 (13) 



B{B^ -> XciK^) 
Combining all of these measurements yields 
B{B'^ —f Charmonium) 



B{B^ -^ Charmonimn) 



1.17 ±0.07 ±0.04 (14) 



Assuming equal partial widths for B'^ — > J/ijjh^ and 
_B^ -^ J/iph'^ for any meson h and using the known ratio 
of the charged to neutral B meson lifetimes ts+ /t^o = 
1.062 ±0.029 |9), we find: 

i?+/o.mi^l^l^^. 1.10 ±0.06 ±0.05 (15) 

We provide the formulae for recomputing our results 
for an arbitrary value of i?+/°, rather than the value of 
unity we have assumed: 

6(B+->X,i?+/°) ^ ii±-^6(S+^X,l)(16) 
6(B" ^ X, i?+/°) = ii±^^6(S° ^ X, 1) (17) 



We also determine the ratio of branching fractions for a 
vector versus scalar light meson accompanying the char- 
monium meson: 



B{B° -^ J/^K*°) 

B{B^^ J/il!K°) 

B{B+ -^ J/tljK*+) 

B{B+ ^J/il)K+) 

BjB^ ^ Xci^*°) 
B{B^ ^ XciKO) 



= 1.49 ±0.10 ±0.08 (18) 
= 1.37 ±0.10 ±0.08 (19) 
= 0.89 ±0.34 ±0.17 (20) 



These three ratios are consistent and yield an average 
value: 



B{B ^ charmonium + vector) 
B{B —I- charmonium + scalar) 



1.40 ±0.07 ±0.06 



(21) 
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TABLE X: The observed number of events in the signal region, estimated background, efficiency, efficiency times secondary 
branching fractions and measured branching fraction for exclusive decays of B mesons involving charmonium. The combinatorial 
background is estimated from a fit to the signal plus sideband region in jties, while the peaking background is estimated with 
Monte Carlo. For the B" — > J/xjjKI mode the inclusive charmonium background is listed in the "Peaking" column and the 
other backgrounds in the "Combinatorial" column. For the branching fractions, the first error is statistical and the second 
systematic. 



Channel 




Nohs 


Combinatorial Bkgr 


Peaking Bkgr 


Efiiciency (%) 


Eflt xS{%) 


Branching fraction/ 10 * 


B'' ^J/^K" 


K'i ^ 7r+^" 


275 


6.1 


± 


2.7 


3.4 ± 


1.1 


33.8 


1.37 


8.5 


± 


0.5 


± 0.6 




7^0 ^ TfO ^« 


77 


12.2 


± 


3.7 


2.3 ± 


0.9 


15.5 


0.29 


9.6 


± 


1.5 


± 0.7 




Kl 


408 


25 


± 


3 


200 ± 


14 


22.3 


1.46 


6.8 


± 


0.8 


± 0.8 




All 


















8.3 


± 


0.4 


± 0.5 


B+ -^ J/i,K+ 




lf35 


8.9 


± 


2.6 


17.1 ± 


2.6 


41.2 


4.86 


10.1 


± 


0.3 


± 0.5 


B^^J/^^TV" 




19 


4.7 


± 


0.9 


0.7 ± 


0.1 


25.8 


3.01 


0.20 


± 


0.06 


± 0.02 


B''^J/i,K*'' 




695 


50.2 


± 


7.8 


50.0 ± 


3.3 


22.6 


1.10 


12.4 


± 


0.5 


± 0.9 


B+ -^ J/i!K'+ 




625 


160.6 


± 


15.9 


87.0 ± 


5.8 


17.9 


1.09 


13.7 


± 


0.9 


± 1.1 


B° -^ ij{2S)K^' 




63 


6.0 


± 


3.3 


1.0 ± 


0.8 


22.0 


0.37 


6.9 


± 


1.1 


± 1.1 


B+ -^ ip{2S)K+ 




247 


27.2 


± 


5.5 


12.5 ± 


2.8 


29.6 


1.46 


6.4 


± 


0.5 


± 0.8 


B -^ XciJ< 




37 


7.2 


± 


2.1 


3.7 ± 


1.3 


19.1 


0.21 


5.4 


± 


1.4 


± 1.1 


B+ ^ XciK+ 




179 


24.2 


± 


4.7 


9.7 ± 


2.7 


26.3 


0.85 


7.5 


± 


0.8 


± 0.8 


B° - Xcii^*" 




52 


13.0 


± 


1.6 


6.4 ± 


5.8 


13.9 


0.30 


4.8 


± 


1.4 


± 0.9 



Finally, the following ratios between the production 
rates for different charmonium states have been deter- 
mined: 



BjB" -^ iIj{2S)K^) 

B{B^ ^ XciK°) 

B{B"^J/i;KO) 

B{B+ -^ tlj{2S)K+) 

B{B+ ^ J/i,K+) 

B{B+ ^ XciK+) 

B{B+ ^J/i>K+) 



0.82 ±0.13 ±0.12 (22) 

0.66 ±0.11 ±0.17 (23) 

0.64 ±0.06 ±0.07 (24) 

0.75 ±0.08 ±0.05 (25) 



XIII. SUMMARY 

We have presented measurements of branching frac- 
tions of B mesons to several two-body final states that 
include a J/ip ,ip{2S) or Xci meson and a K'^ ^ K'^ ^ K* 
or TT*^. Our results are in good agreement with previ- 
ous measurements [|l9| and have superior precision, both 
in terms of individual branching fractions and their ra- 
tios. In addition, based on isospin invariance, we find the 
ratio of charged to neutral B meson production on the 
T(4S') resonance to be compatible with unity within two 
standard deviations, and also compatible with the mea- 
surement reported by CLEO |26| . Our central value and 
CLEO's are both higher than one, with the difference in 
our case larger than one standard deviation. 
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